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Abstract. Properties of inertial modes of Jupiter are investigated for an n = 1 polytropic description of the planet 
interior. We use the anelastic approximation to overcome the usual handicap of a severe spherical harmonics 
truncation. A powerful iterative solver then allows us to compute the frequencies of the most promising low-order 
modes using as many spherical harmonics as necessary. The induced 0(1%) errors of our model are now within 
observational limits. A plausible seismological model of Jupiter might thus be at hand provided the use of a more 
realistic description of the planet interior. 
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1. Introduction 

Detecting global Jovian oscillations is still a major chal- 
lenge since, even with the most recent detection tech- 
niques, no acoustic or gravity modes have been con- 



vinci n gly id entified (Schmi der e t al. 1991; Mosser et al. 
1993| , ^OOOj Cacciani et al. |200lD . Furthermore, the colli- 



sion of the Shoemaker-Levy9 fragments was unfortunately 
not energetic enough as to al low an y positive detection of 
primary waves (Walter et al. 1996 Mosser et al. 1996 ). 
Upcoming observational projects, such as JOVIS 



(Baglin & Mosser 1999), open a new window on the planet 



dynamics with the possibility of detecting long-period 
mode oscillations. Because of Jupiter's rapid rotation, in- 
ertial modes would be good candidates. These have first 
been investigated by Lee & Saio (1990) and Lee et al 1992 



(hereafter LSVH) using different models for the planet. 
Their approach was nonetheless limited by the use of no 
more than three spherical harmonics in the angular de- 
scription of the solutions (see § [2.l| ). We have today the 
numerical tools to tackle the same problem by including 
as many terms as necessary in each of the spherical har- 
monic expansions (Rieutord & Valdettaro 1997; Dintrans 
et al. |1999D. 



The plan of the paper is as follows: in §H we intro- 
duce our model, its specificity and the equations we solve 
for. In we discuss mode identification methods in pres- 
ence of rotation and show how the inclusion of diffusion 
allows us to discriminate the inertial eigenmodes in terms 
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of least/high damped ones. Our results are then presented 
in §^ before concluding in 



2. Our model 

2.1. Specificity of our study 

The properties of inertial mode oscillations are commonly 
investigated by first expanding each perturbation in a serie 
of spherical harmonics. However, because of the Coriolis 
couplings between the harmonics (^,£± 1), one is faced 
with an infinite system of coupled ordinary differential 
equations. Such a system is challenging in its coding and 
requires prohibitive memory. LSVH solved these two as- 
pects by considering at most the first three terms in each of 
the spherical harmonic expansions. They nevertheless rec- 
ognized the necessity of calculations with larger numbers 
of spherical harmonics before detailed comparisons could 
be made between the theoretically calculated frequencies 
and the yet to be measured oscillations of Jupiter. Indeed, 
for the higher overtones, truncation errors as high as 30% 
were reached when two harmonics instead of three were 
used. 

In our case, the memory requirement is immediately 
dealt with by adopting the anelastic approximation (see 
also Dintrans & Rieutord 200C). The idea is to filter out 
the acoustic waves in the infinite system while keeping 
the density variations across the planet. This is further 
justifiable since such high-frequency waves are irrelevant 
to the dynamics of the inertial modes. Remains the cod- 
ing aspect of the problem which is taken care of by a 
preprocessing Perl program capable of managing as many 
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spherical harmonics as necessary. Errors in our approach 
are thus mainly induced by the use of the anelastic ap- 
proximation and are estimated to be 0(1%) (see Dintrans 
& Rieutord |00l|). 



2.2. The oscillation equations under the anelastic 
approximation 

Assuming the time-dependence of the eigenmodes to be 
expiat (where a is the angular frequency in units of 
rad/s), the anelastic equations for the linear perturbations 
in the co-rotating frame are given by 



iau + 2n X M = -V(P7po), 
div{pau) = 0, 



(1) 



where u is the velocity, P' the eulerian pressure perturba- 
tion, fl the angular frequency (217 x u being the Coriolis 
force) and po denotes the equilibrium density of our in- 
terior model of Jupi ter wh ich we describe by the n = I 
polytrope (Hubbard |l984 ). As boundary conditions, we 
impose the regularity of the velocity at the centre while 
the anelastic approximation implies that the radial com- 
ponent of th e velo city vanishes at the surface (see Dintrans 
& Rieutord ^00l|) . 

Using the spherical coordinate system (r,0,ip), we 
expand t he ve locity on spherical harmonics Y™{9,(p) 
(Rieutord 1991 ). One arrives at an infinite set of coupled 
radial equations to be solved as a generalized eigenvalue 
problem of the form 



1±, 



(2) 



where X — ia is the complex eigenvalue associated with 
eigenvector ip„i± . The m^-modes corresponds to even 
solutions {£ = |m|,|TO| -I- 2, |m| -I- 4,...) whereas m~- 
modes are odd ones {£ = |m| -I- 1, \m\ + 3, \m\ -I- 5, . . .). 
We note that even (odd) modes define solutions that are 
symmetric (antisymmetric) with respect to the equatorial 
plane. Finally, this eigenvalue problem is discretized on 
the Gauss-Lobatto grid associated with Chebyshev poly- 
nomials yielding to matrices A4a and A^s of order about 
L X (A'^ -|- 1) when L spherical harmonics and A^ polyno- 
mials are considered. 

3. Mode identification in presence of rotation 

In this section we discuss the notion of mode identification 
when adding rotation to the p- and g-modes as compared 
to the case of inertial modes. 



3.1. p- and g-modes 

One recalls that two cases need to be considered depending 
on the non-rotating frequency shift aki {k being the radial 
order) with respect to the Coriolis frequency 257: 
i) cTke ^ 2r2 {p- and low-order g-modes): the weak Coriolis 
couplings can be described by a perturbative theory (i.e. 



akim — ^ki + mCke^) making the {k,£,m) identification 
possible. 

ii) crfe£<,251 (high-order g-modes): the strong Coriolis cou- 
plings prevent any ^-identification of modes. A pulsation 
mode cannot henceforth be identified by a set (£, m) since 
the description of its angular dependence in terms of a sin- 
gle spherical harmonic Yl^{6, Lp) is not possible (Dintrans 
et al. 1999). 

3.2. Inertial modes 

Only the aki < 257 regime exists, i.e. the spectrum is 
bounded by the Coriolis frequency. The Coriolis couplings 
are dominant over the entire spectrum making the £- 
identification meaningless. Although they recognized this 
fact, LSVH chose "for simplicity" to attribute a given 
^-value to their inertial eigenmodes; i.e. they arranged 
their low-order solutions in terms of £ = 1,2,3,4 for 
m — —1, —2, and m = —3 with eigenvalues sorted 
in decreasing values. This arrangement is in fact unap- 
propriate, as shown for example with the eigenmodes 
of the incompressible rotating sphere (Bryan 1889): for 
m = 0"*", the first four low-order frequencies are a/2Q = 



[0.654; 0.830; 0.469; 0.677' 



(see also Greenspan 196E) 



with no possible ^-attribution 



3.3. Identification from the (I'll, a)-diagram 

Low-degree and low-order p- or g-modes are commonly 
identified in the computed spectrum by plotting the {£, a)- 
diagram. Unfortunately, such a tool cannot be used for 
inertial modes since £ is not a suitable quantum number 
anymore. 

However, one would expect that the most promising 
observational inertial modes are the ones with the weak- 
est Coriohs couplings. Such modes, hereafter referred to as 
"low-order" modes, only involve a small family of spher- 
ical harmonics and are for this reason associated to the 
smallest resolutions L. Therefore the challenge is to ex- 
tract them among the entire computed modes defining 
the spectrum, thus the question: "At a given m and fre- 
quency shift cr, what is the inertial mode which involves 
the least of spherical harmonics, that is, the one whom 
the ^-couplings are the weakest ?". As we will prove in 
§4.2 and in more details in a forthcoming paper, one way 
to extract the modes in question consists in plotting a 
(|x|, (T)-diagram, x being the (negative) damping rate of 
the mode when adding a slight diffusion. 

4. Results 

Hereafter dimensionless angular eigenfrequencies a and 
damping rates T are given in units of 2f2, that is ct = 2Qa 
and Z— 2^lz. For Jupiter, it means that oscillation fre- 
quencies are i^(mHz) = a/2Tr ~ 0.056ct whereas periods 
are T(hr) = 2Tr/a ~ 4.958/5. 
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Fig. 1. The (|t|, cr)-diagram for the m — 0^ eigenvalues, 
where x and a are the dimensionless damping rate and 
frequency, respectively. The thick dots to the left denote 
the first three low-order modes. The diamond, square and 
triangle modes, chosen in order of increasing ^-couplings 
(arrow), are used when discussing truncation effects in 



4.1. Low-order modes for m = — 1+, O"*", l"*" 

We add a small viscosity term to the momentum equation 
in order to organize the eigenvalues in terms of least/high 
damped pulsation modes. By plotting them in the (|x|, cr)- 
diagram, we are able to identify the low-order modes 
which position themselves in the left region of the plot. 
An example is shown in Fig. for the m = eigenval- 
ues where few low-order modes are depicted by thick dots. 
This figure was obtained using the QZ algorithm with the 
resolution L = N = 60 and an Ekman number E — 10^^ 
{E = v/2Q.R^ , where v is the kinematic viscosity and R 
the Jupiter radius). 

Once the best low-order modes isolated, remains the 
important task/step of accurate computation of the cor- 
responding eigenfrequencies in the astrophysical relevant 
limit of zero viscosity. The QZ algorithm is not the most 
suitable for this job because the matrix band structure 



Table 1. Dimensionless frequencies a of the first five low- 
order modes for m = — 1+, 0+, 1+. The associated oscilla- 
tion periods (in hours) are given in parenthesis. 



0.656 (7.558) 
0.201 (24.67) 
0.799 (6.205) 
0.491 (10.10) 
0.327 (15.16) 



0+ 

0.679 (7.302) 
0.828 (5.988) 
0.507 (9.779) 
0.403 (12.30) 
0.893 (5.552) 



0.706 (7.023) 
0.863 (5.745) 
0.719 (6.896) 
0.921 (5.383) 
0.343 (14.45) 



Fig. 2. Normalized angular power spectra of the radial 
component Ur of the velocity for the three modes labelled 
by a diamond, a square and a triangle in Fig. ^ Note the 
increase of the ^-couplings with the mode order. 



of M^A A4g is not taken advantage of; i.e. Fig. |^ re- 
quired 1.8 Gbytes of memory and 16 hours of computation 
on Origin2000. 

We use an Arnoldi-Chebyshev iterative solver to ex- 
tract the interesting eigenvalues lost in the dense spec- 
trum at i5 = 0. Taking advantage of the matrix band 
structure, it leads to substantial gain of both memory 
and computing time. Table |l| shows the computed fre- 
quencies of the first five low-order modes for the m — 0^ 
and m = ±1+ cases. As can be noted from the associated 
periods given in parenthesis (in hours), low-order inertial 
modes - which must always exceed T^ot/'^ ^ 4.958 hr - can 
be as high as one day. It is thus clear that these very long 
period modes require continuous observation using either 
a ground-based network of dedicated telescopes or, yet 
better, a devoted satellite such as the mentioned JOVIS 
project. 

4.2. Truncation effects 

In order to test the effect of the harmonic truncation, we 
compute, in the limit E — the ^-spectra of the three 
modes denoted by a diamond, a square and a triangle in 
Fig. |l|. That is, for each £, the maximum value over the 
Chebyshev coefficients of Ur is plotted. 

Fig. 1^ shows that only few spherical harmonics are nec- 
essary to describe the lowest-order eigenmode (diamond) 
whereas many more harmonics are required to resolve 
the square and triangle ones. To reach a relative error 
of 0(10"^), we indeed checked that the diamond mode 
only requires L — A spherical harmonics while the square 
and triangle modes necessitate at least L = 14 and L = 22 
harmonics, respectively. That is, eigenvalues are extremely 
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sensitive to L-values before the sharp drop in the curves 
of Fig. I 

Combining our findings from Figs. ^ and ^, the Hnk 
between the damping rate of an inertial mode and the 
strength of its ^-couphngs is clearly demonstrated: the 
higher the damping rates, the higher the ^-couplings, thus 
the higher the needed resolution L. It well justifies the 
use of the (|t|, CT)-diagram to identify the low-order iner- 
tial modes and to remedy to the fact that £ is no longer a 
suitable quantum number. 

5. Conclusion 

We investigated the properties of the inertial modes of 
Jupiter for an n = 1 polytropic model of the planet. Our 
calculations have been performed with as many spheri- 
cal harmonics as necessary to resolve eigenmodes in the 
anelastic approximation. 

We showed that adding a slight viscosity to the equa- 
tions allows us to identify what we referred to as "low- 
order" modes, that is inertial modes which only involve a 
few spherical harmonics and are the most promising for 
observations. Once these low-order modes identified, we 
compute them in the relevant adiabatic limit of zero vis- 
cosity. We also showed that, while the lowest-order mode 
is fully described by four spherical harmonics, the follow- 
ing low-order modes require many more harmonics (up to 
twenty) as a consequence of the succeeding Coriolis cou- 
plings. 

An avenue for future work consists on taking into ac- 
count density discontinuities in the planet's interior (such 
as the so-called PPT or the plasma phase transition; 
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